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ABSTRACT 

We present an analysis of F-band radial surface brightness profiles for spiral galax- 
ies from the field and cluster environments using Hubble Space Telescope/ Advanced 
Camera for Surveys imaging and data from the Space Telescope A901/2 Galaxy Evo- 
lution Survey (STAGES). We use a large sample of ~ 330 face-on to intermediately 
inclined spiral galaxies and assess the effect of the galaxy environment on the az- 
imuthally averaged radial surface brightness /i profiles for each galaxy in the outer 
stellar disc (24 < /i < 26.5 mag arcsec -2 ). For galaxies with a purely exponential 
outer disc (~ 50 per cent), we determine the significance of an environmental depen- 
dence on the outer disc scalelength h OU f For galaxies with a broken exponential in 
their outer disc, either down-bending (truncation, ~ 10 per cent) or up-bending (anti- 
truncation, ~ 40 per cent), we measure the strength T (outer-to-inner scalelength 
ratio, log 10 h out /hi n ) of the /i breaks and determine the significance of an environmen- 
tal dependence on break strength T. Surprisingly, we find no evidence to suggest any 
such environmental dependence on either outer disc scalelength h ou t or break strength 
T, implying that the galaxy environment is not affecting the stellar distribution in the 
outer stellar disc. We also find that for galaxies with small effective radii (r e < 3kpc) 
there is a lack of outer disc truncations in both the field and cluster environments. 
Our results suggest that the stellar distribution in the outer disc of spiral galaxies is 
not significantly affected by the galaxy environment. 
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1 INTRODUCTION 

The structure of the outer regions of galactic discs is im- 
portant to our understanding of the formation and evo- 
lution of spiral galaxies. These faint, outer regions are 
more easily affected by interactions with other galaxies, and 
therefore their structural characteristics must be closely re- 
lated to their evolutionary history. Certain physical pro- 



cesses inherent to galaxy evolution and dependent on 
the galaxy environment, for example, ram pressure strip- 
ping of the interstellar med i um, mergers, and harassment 
(e.g.. iGunn fc Gottl If 9721 : lOkamoto fc Nagashimal 12004 
M oore et al.lll996l ) could have an effect on the galactic disc 
and therefore the light distribution (surface brightness \x 
profile) of a galaxy. Therefore, the environmental depen- 
dence of the shape of spiral galaxy radial surface brightness 
profiles in the outer galactic disc will provide evidence for 
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the physical processes of galaxy evolution occurring in dif- 
ferent environments. 

Since t he work of Patterson! (|l940l ). Ide Vaucouleursl 

l|l959l ). and iFreemar] l|l970l ) we have known that the light 
profiles of spiral galaxies are comprised of two main com- 
ponents: an inner component dominated by a bulge, and 
an outer disc that follows a simple exponential decline with 
minor deviations related to the spiral arms. However, this 
classical picture does not hold for most spiral galaxies in 
the Universe and has been shown to f ail at the faint sur- 
face brightness of t he outer stellar disc dvan der Kruitll 19791 ; 
iPohlen et all |2002|) . Since Ivan der Kruitl l| 19791 ) we have 
known that the exponential decline in the outer disc does not 
extend out to the last measured point but can be truncated 
(sharply cut off) a fter several scalelengths. More recently 
IPohlen et all (|2002l ) have shown that the exponential de- 
cline does not cut off completely at the truncation. They find 
that most profiles are actually best described by a two slope 
model (broken exponential), characterised by an inner and 
outer exponential scalelength separated by a relatively well 
defined break radius rb r k- Many studies have now reported 
(mainly using surface photometry) the existence of broken 
exponential ste llar discs (truncations) in spiral galaxies in 
both the local (IPohlen et all |2002|; IPohlen fe Truiilloll2006l ; 
IPohlen et ai]|2007l ; Ivan der Kruitl 120071; iBakos et alj 120081) 



may account for a sharp truncation (cut off), but does not 
explain observations of exten ded outer exponent ial compo- 
nents in some galaxies (e.g. IPohlen et al.l |2002|). An other 



theory has been proposed by Dcbattist a et al. (2006) who 



and distant z < 1 Univers e l|Pereal2004l ; iTruiillo fe Pohlenl 
120051 ; I Azzollini et al.ll200i ). 

As a result of these studies it has been shown that 
disc galaxies can be classified into three broad types (Type 
I, II, and III) according to break featu res in their radial 
stella r surface brightness profiles (see e.g. IPohlen fe Truiilicl 
2006). Type I (no break), the galaxy has a simple expo- 
nential profile extending ou t to several scalelengths (e.g. 
iBland-Hawthorn et al.ll2005l ). Type II (down-bending break, 
truncation), the exponential is broken with a shallow in- 
ner and steeper outer exponential re gion separated by a 
relatively well defi ned break radius l|van der Kruitl fl979l ; 
IPohlen et al. 1 120021 ). Type III (up-bend ing break, antitrun - 
cation), a class recently discovered by lErwin et all l|2005l ) 
which have a broken exponential with the opposite be- 
haviour to a Type II profile. In each case the classification 
refers to the outer, disc component of the galaxy radial sur- 
face brightness profile and does not consider the inner vary- 
ing bulge component. In some cases the inner bulge com- 
ponent may be near exponential in nature, however, this 
classification scheme only considers the disc component and 
is fairly insensitive to the nature of the inner (bulge) profile. 

Measurements independent of surface photometry 
(from resolved star counts) are also available on 
nearby galaxies for each of the three profile types. 
IBland-Hawthorn et all (|2005l ) find that NGC300 has a sim- 
ple expo nential profile extend ing out to ~ 10 scalelengths 
(Type D. lFerguson et all (|2007h argue that M33 is best de- 
scribed as a br oken exponential w ith a down-bending break 
(Type II), and llbata etail (|2005l ) report that M31 could be 
described as having an antitruncated disc (Type III). 

The physical processes that cause the different types 
of profile are not well understood. Some models suggest 
that Type II profiles (truncations) could be due to the 
effect of a star formation threshold in disc column den- 
sity towards the edge of the disc (e.g. | Kennicufrd Il989l ; 
lElmegreen fe Parravand 1 1994 ISchave! I2004T ). This theory 



find down-bending breaks in simulated disc profiles solely 
from collisionless N-body simulations. They suggest Type II 
galaxies are the consequence of a resonance phenomenon and 
a redistribution of angular momentum that leads to an in- 
creased central density and surface brightness. Many models 
now incorporate both these ideas and suggest that the inner 
disc forms as a consequence of the star formation threshold 
while the outer disc forms by the outward migration of stars 
from the inner disc to regions beyond the star formation 
threshold. This migratio n could be due to r esonant scat- 
terin g with spiral arms (IRoskar et al. 2008a: IRoskar et al.l 
12008b ) or clump disruptions (|Bournaud et al.ll2007l) . Alter - 
natively, another model is suggested by iFovle et all (|2008l ) 
who use simulations to show that by starting from a single 
exponential disc, the inner disc forms as the bulge draws 
mass from the inner regions. As a consequence the inner 
disc profile becomes shallower while the outer region stays 
almost unaffected. Less work h as been carried out o n the 
origin of Type III galaxies, but lYounger et ail (|2007l ) have 
shown that recent minor mergers could produce up-bending 
stellar profiles in the remnant galaxy. 

However, all these models rely on the break in the sur- 
face brightness profile having an analogous break i n the 
stellar surface mass density profile. IBakos et al l (|200Sl ) used 
colour profiles to calculate stellar mass surface density pro- 
files for a sample of Typ e II an d Type III galaxies identi- 
fied by IPohlen fe Truiilicl (|2006l ). Interestingly, for Type II 
galaxies they find the stellar mass surface density profiles 
are almost purely exponential. Therefore, the break in the 
surface brightness profile is not necessarily related to a break 
in the stellar mass surface density profile. Consequently, the 
surface brightness break may not be due to the distribution 
of stellar mass but could be due to a radial change in the 
stellar population. For Type III galaxies however, they do 
find an analogous break in the stellar surface mass density 
profile indicating these breaks could be due to the stellar 
mass distribution. 

Investigating the frequency of profile types in differ- 
ent galaxy environments will provide evidence for their 
origin and the effect of the environment on the outer 
stellar disc. Presently, there have been few systematic 
searches for stellar disc truncations in jspiral galaxies 



in either the local dPohlen fe Truiilld I2006T) or distan t 



Universe (Truiillo fe Pohlenl 120051 ; I Azzollini et al.l [2008) 



IPohlen fe Truiilicl ((2006) use a sample of ~ 90 face-on to in- 
termediate inclined nearby late-type (Sb-Sdm) s piral galax- 
ies fr om the Sloan Digital Sky Survey (SDSS; lYork et al.l 
2000) and find that approximately 10 per cent are Type I, 
60 per cent are Type II and 30 per cent are Type III. They 
also report that the shape of the profiles correlate with Hub- 
ble type. In their sample down-bending breaks are more fre- 
quent in later Hubble types while the fraction of up-bending 
breaks rises towards earlier types. They also find no rela- 
tion between the galaxy environment, as determined by the 
number of nearest neighbours, and the shape of the surface 
brightness profile. However, they do not reach very dense 
environments and low number statistic s did not allow for 
major conclusions. I Azzollini et all (|2008l ) recently conducted 
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the largest systematic search for stellar disc truncations yet 
undertaken at intermediate redshift (0.1 < z < 1.1) us- 
ing the Great Observa tories Origins Deep Survey (GOODS; 
iGiavalisco et al.ll2004T ) south field. They use a sample of 505 
galaxies and obtain the frequency of profile types in differ- 
ent redshift ranges. They find that the frequency of profile 
types (Type 1:11:111) is 25:59:15 per cent for 0.1 < z < 0.5, 
and does not vary significantly with redshift out to z ~ 1.1. 

The aim of this paper is to undertake a systematic 
search for broken exponentials in the field and cluster en- 
vironment using the Spa ce Telescope A901 /2 Galaxy Evo- 
lution Survey (STAGES: [Gray et all 120091 ') and to investi- 
gate whether the galaxy environment has an effect on the 
stellar distribution in the outer stellar disc. We investigate 
whether the type and strength of radial /x profile breaks in 
the outer stellar discs of spiral galaxies is dependent on the 
galaxy environment. This work builds on previous studies 
by using larger and more statistically viable field and clus- 
ter samples and by reaching higher density environments. 
However, it is important to note that STAGES only cov- 
ers an intermediate density environment (projected galaxy 
numb er density up to ~ 1600 galMpc -3 . Irleiderman et al.l 
2009) and not a high de nsity environment (e .g. the COMA 
cluster, ~ 10 4 galMpc~ 3 . lThe fc Whitl 19861 ). We also wish 
to point out that in this work we use a sli ghtly different pro- 
file cl assific ation scheme to that u sed bv lPohlen fc Truiillol 
(2006) and lAzzollini et al l (|2008h . see Section 4.3. There- 
fore, direct comparisons of the results of this work to the 
previous works mentioned above cannot be made. 

The structure of this paper is as follows: in Section[2]we 
give a brief description of the STAGES dataset relevant to 
this work and outline our sample selection in Section \2. II In 
Section [3] we describe the method used to obtain our radial 
surface brightness profiles from the STAGES Hubble Space 
Telescope (HST)/ Advanced Camera for Surveys (ACS) V- 
band imaging and explain our profile classification scheme 
in Section [4] We present our results in Section [5] and finally 
draw our conclusions in Section [B] Throughout this paper, 
we adopt a cosmology of Ho = 70kms~ 1 Mpc~ 1 , = 0.7, 
and £l m — 0.3, and use AB magnitudes unless stated other- 
wise. 



2 DESCRIPTION OF THE DATA 

T his work is entirely based on the STAGES data published 
bv lGrav et"afl (12009T ). STAGES is a multiwavelength survey 
that covers a wide range of galaxy environments. A com- 
plex multicluster system at z ~ 0.167 has been the sub- 
ject of V-band (F606W) HST/ ACS imaging covering the 
full 0.5° x 0.5° (~ 5 x 5Mpc 2 ) span of the multicluster sys- 
tem. The ACS imaging is complemented by photometric red- 
shifts and observed- /rest-frame spectral energy distributions 
(SEDs) from the 17-ban d COMBO-17 photometric redshift 
survey (|Wolf et alj|2003r ). Extensive multiwavelength obser- 
vations using Spitzer, Galaxy Evolution Explorer (GALEX), 
2 degree Field (2dF), XMM-Newton, and the Giant Metre- 
wave Radio Teles cope (GMRT) have also been carried out. 
iGrav et al.l (|2009l) have performed Sersic profile fitting using 
the GALFIT code l|Peng et alj [2002) on all HST/ ACS images 



and conducted simulations to quantify the completeness of 
the survey, all of which are publicly availablfl 

The COMBO-17 observations used in the STAGES 
master catalogue were obtained with the Wide Field Imager 
(WFI) at the Max Planck Gesellschaft/European Southern 
Observatory (ES Q) 2.2-m telescope on La Silla, Chile (see 
IWolf et all [2003! for further details). COMBO-17 used five 
broad-band niters UBVRI and 12 medium-band niters cov- 
ering wavelengths from 350-930 nm to define detailed opti- 
cal SEDs for objects with R ^ 24, with 7? being the total 
Vega i?-band magnitude. Generally, photometric redshifts 
from COMBO-17 are accurate to 1 per cent in 5z/(l + z) at 
R < 21 which has been spectroscopically confirmed. Photo- 
z quality degrades for progressively fainter galaxies reaching 
accuracies of 2 per cent for gala xies with R ~ 22 and 1 per 
cent for galaxies with R > 24 (|Wolf et al.ll2004 jpoj ). To 
date, all galaxy evolution studies on the COMBO-17 data 
that use photo- z defined galaxy samples restrict themselves 
to galaxies with R < 24 to ensure only reliable redshifts 
are used. Stellar mass estimates derived from SED fitting 
the 17-b and photometry a r e also available f or COMBO-17 
galaxies jBorch et al.ll2006l : IGrav et alj|2009h . 

The STAGES morphological catalogue (Gray et al. in 
prep.) contains 5090 galaxies in STAGES with reliable Hub- 
ble type morphologies. All galaxies with R < 23.5 and 
Zphot < 0.4 were visually classified by seven members of 
the STAGES team into the Hubble types (E, SO, Sa, Sb, 
Sc, Sd, Irr) and their intermediate classes. SOs were defined 
to be disc galaxies with a visible bulge but no spiral arms 
(smooth disc). Spiral galaxies have visible spiral arms in the 
disc and the classification (Sa-Sd) represents a decreasing 
sequence in the bulge-to-disc (B/D) ratio. Weighted aver- 
age estimates of the Hubble types, ignoring bars and degrees 
of asymmetry were generated. In this paper, we only con- 
sider visually classified spiral galaxies (Sa, Sb, Sc, Sd) and 
intermediate spiral classes are grouped to the earlier Hubble 
type (e.g. Sab are considered to be Sa). 

2.1 Sample selection 

IGrav et all l|2009f ) suggest a cluster sample for STAGES de- 
fined solely from photometric redshifts. The photo- z distri- 
bution of cluster galaxies was assumed to follow a Gaus- 
sian, while the field distribution was assumed to be con- 
sistent with the average galaxy counts N(z, R) outside the 
cluster and to vary smoothly with redshift and magni- 
tude. The cluster sample is defined by a redshift interval 
Zphot = [0.17— Az, 0.17 + Az], where the half- width Az was 
allowed to vary with 7?-magnitude. A narrow redshift range 
is adopted for bright i?-magnitudes due to the high preci- 
sion of COMBO-17 photometric redshifts, but the interval 
increases in width towards fainter ^-magnitudes to accom- 
modate for the increase in photo-z error. The completeness 
and contamination of the cluster sample at all magnitude 
points was calculate d using the coun ts of the smooth mod- 
els (see Fig. 14 from lGrav et alj|2009h and the photo-z half- 
width compromised so the completeness was > 90 per cent 
at any magnitude. Contamination (< 25 per cent) is defined 
to be the fraction of field galaxies in the cluster sample at 

1 http : / / www.nottingham. ac . uk/astronomy /stages 
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a given magnitude (and not below it). The half- width as a 
function of magnitude R is 



Az(R) = v/0.015 2 + 0.0096525 2 (1 + 10°- 6 <^ 



-20.5) ) 



(1) 



This equation defines a photo-z half-width that is lim- 
ited to 0.015 at bright ^-magnitudes but increases as a 
constant multiple of the estimated photo-z error at the 
faint end. The completeness of this sele ction converges to 
nearly 100 per cent for bright galaxies; see lGrav et al. (l2009h 



0.22,0.30] 
Wolf et all 



for fu rther details. The catalogue published by iGrav et al] 
(2009) contains a number of flags (combo_flag, phot_flag, 
stages_flag) that allows the selection of various galaxy sam- 
ples. In this study, we use a similar sample selection to that 
used bv lMaltbv etail (|2010h . 

For our cluster sample we use the above cluster defini- 
tion (combo_flag ^ 4), we also only use galaxies with reli- 
able photometry (i.e. those with phot_flag < 8), and those 
defined as extended HST sources in STAGES (stages_flag 
^ 3). We also limit our sample by stellar mass, cutting at 
logAf, /Mq > 9, and only select galaxies with visually clas- 
sified Hubble-type morphologies in the range Sa-Sdm us- 
ing the STAGES morphological catalogue (Gray et al. in 
preparation) . This cluster spiral galaxy sample contains 383 
galaxies. 

For our field sample we use COMBO-17 defined galaxies 
(combo_flag ^ 3) and apply a redshift selection that avoids 
the cluster. We include a lower redshift interval at z = 
[0.05, 0.14] and an upper redshift interval at z = 
based on a similar sample selection used by 
(2009). We also only use galaxies with reliable photome- 
try (i.e. those with phot_flag < 8), and those defined as ex- 
tended HST sources in STAGES (stages.flag > 3). We also 
limit our sample by stellar mass, cutting at log M,/Mg > 9, 
and only select galaxies with visually classified Hubble-type 
morphologies in the range Sa-Sdm using the STAGES mor- 
phological catalogue (Gray et al. in preparation). This field 
spiral galaxy sample contains 31 8 galaxies. 

The catalogue published bv lGrav et al] (2009) contains 
two sets of derived values for properties such as magnitude 
and stellar mass, one based on the photo-z estimate and 
another assuming the galaxy is located at the cluster redshift 
of z = 0.167. This prevents the propagation of photo-z errors 
into physical values. Here we use the fixed redshift set of 
values for the cluster sample, but the original estimates for 
our field comparison sample. 

The field and cluster samples are also restricted by 
galaxy inclination i to select face-on to intermediately in- 
clined spiral galaxies by using the axis ratio q o f the galaxy 
as de termined by the STAGES galfit models (|Grav et al] 
2009). The minor-to-major axis ratio q (q — b/a = 1 — e, 
where a and b are the semi-major and semi-minor axes 
respectively and e is the ellipticity) is restricted to corre- 
spond to an inclination i of less than 60 degrees (q > 0.5 
or e < 0.5). This cut is necessary to avoid the influence of 
dust on our surface brightness profiles and allows reliable in- 
formation on features like bars, rings, and spiral structure. 
The distribution of axis ratio q for our field and cluster spiral 
samples showing the axis ratio q cut is shown in Fig. [T] This 
cut removes ~ 54 per cent of field spirals and ~ 52 per cent 
of cluster spirals from our samples. The final cluster sample 
contains 182 spiral galaxies reaching down to R ~ 22 and 
has a photo-z range of z p hot = [0.143, 0.196]. The final field 
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Figure 1. The distribution of minor- to-major axis ratio q for 
our field (light grey) and cluster (dark grey) spiral samples. The 
q cut (q > 0.5, represented by a black dashed line) ensures the 
selection of face-on to intermediately inclined spiral galaxies (i < 
60 degrees). Errors in q are < 3 per cent. Respective sample sizes 
are shown in the legend. 
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Figure 2. The photometric redshift estimate versus total R- 
band magnitude for the final field (grey points) and cluster (black 
points) spiral galaxy samples. The field sample reaches R ~ 23 
and the cluster sample reaches R ~ 22. Respective sample sizes 
are shown in the legend. 



sample contains 145 spiral galaxies reaching down to R ~ 23 
(see Fig. [2]). The morphological mix (Sa, Sb, Sc, Sd) of the 
final field and cluster spiral samples are shown in Fig. [3] 

STAGES is > 90 per cent complete for R < 23.5 
l|Grav et al] |2009| ) as is the case for our field and cluster 
galaxy samples (see Fig. [2]). However. IWolf et al] (|2009T ) es- 
timate that at stellar masses below logM*/M0 < 9.5, the 
field sample could have an additional 20 per cent incomplete- 
ness based on previous COMBO-17 experience. Therefore, 
our field sample is essentially > 70 per cent complete. We 
expect the completeness of the cluster sample to be > 90 
per cent and the contamination of the cluster sample by 
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Figure 3. The distribution of Hubble-type morphologies for the 
final field (light grey) and cluster (dark grey) spiral galaxy sam- 
ples. Respective sample sizes are shown in the legend. A signif- 
icant excess of earlier types (Sa) is observed in the cluster envi- 
ronment. 



Table 1. Properties of the field and cluster spiral samples. 



Property 


Field 


Cluster 




145 


182 


Completeness 


> 70 per cent 


> 90 per cent 


Contamination 




< 25 per cent 


■Rmcai! 


20.14 


19.32 


^B(min) 


-15.1 


-16.1 


^B(max) 


-21.5 


-21.8 


z p hot, me an 


0.227 


0.171 


^•phot.min 


0.055 


0.143 


2p hot, max 


0.299 


0.196 



the field to be < 25 per cent, based on the fl -magnitude 
the r espective samples reach (see Fig. 14 from iGrav et al.l 
2009). Selecting spiral galaxies from the STAGES morpho- 
logical catalogue (Gray et al. in prep) i ntroduces no furthe r 
incompleteness to the galaxy samples (|Maltbv et al.l l2010). 
The properties of the final field and cluster spiral samples 
are shown in Table. [T] 



3 PROFILE FITTING 



We use the IRAF task ellipse (STSDAS package - version 
2.12.2) in order to obtain azimuthally averaged radial sur- 
face brightness profiles from the STAGES HST/ACS V- 
band imaging. Profile fitting is carried out for each spiral 
galaxy in our field and cluster samples. The ACS images 
used include the sky background and the necessary sky sub- 
traction is performed after the profile fitting (see Section 
13. 2[) . For further deta i ls on the fitting method used by el- 
lipse see I Jedrzeiewskil (| 19871 ). 

We run ellipse using bad pixel masks that remove 
flagged pixels from the isophotal fit. This is necessary in 
order to remove sources of contamination such as back- 
ground/companion galaxies and foreground stars (every- 



thing not associated with t he galaxy itself), see Fig. [4]for an 
example. IGrav et al.l (|2009l ) generated these bad pixel masks 
for each STAGES galaxy using the data pipeline GALAPAGOS 
(Galaxy Analysis over Large Areas: Parameter Assessment 
by GALFlTting Objects from SExtractor; Barden et al., in 
prep). GALAPAGOS performs an extraction of source galax- 
ies from the STA GES HST/ACS F -band imaging and uses 
the galfit code l|Peng et al.l l2002h to fit ISersid l|l968l ) ra- 
dial surface brightness models to each galaxy image. The 
bad pixel masks are generated for each STAGES galaxy in 
order to remove sources of contamination from these surface 
brightness models. However, occasionally in crowded regions 
GALFIT performs multiobject fitting and therefore the com- 
panion galaxies in these cases are not removed by the bad 
pixel mask as they are too close to the subject galaxy. In 
these cases (27 field and 33 cluster galaxies, ~ 18 per cent) 
the companion galaxies are removed from the ACS image by 
subtraction of their galfit surface brightness model. The 
residuals of the companion galaxies are not expected to have 
any significant affect on the azimuthally averaged radial sur- 
face brightness profile for the subject galaxy. Isophotal fit- 
ting is only attempted when the fraction of flagged (masked) 
data points in the isophote is less than 0.5 and fitting is ter- 
minated if this condition is broken. 

For all our ellipse isophotal fits the galaxy centre is 
fixed (all isophotes have a common centre) using the cen- 
tre of the galaxy determined from the galfit Sersic model 
|Grav et al.ll2009r ). We also use logarithmic radial sampling 
with steps of 0.03 dex (0.07 geometric sampling in ellipse 
terminology) and start from an initial semi-major axis of 10 
pixels. 

For each galaxy in our sample we fit two different sets 
of ellipses to the galaxy ACS image. The first is a free- 
parameter fit (fixed centre, free ellipticity e and position an- 
gle PA) and tends to follow morphological features like bars 
and spiral arms, and therefore is not ideal for the charac- 
terisation of the underlying outer stellar disc studied in this 
paper. Therefore, we use a fixed-parameter fit (fixed centre, 
e, and PA) in order to produce our final surface brightness 
profiles. 

The initial free-parameter fit (fixed centre, free e and 
PA) is used to determine the ellipticity and position angle 
of the outer disc component. For each sample galaxy the 
ACS image was inspected with the contrast adjusted in order 
to get an estimate of the semi-major axis for the end of 
the stellar disc <2di sc Um (where the galaxy surface brightness 
enters the background noise). This 'stellar disc limit' eidisciim 
corresponds to the outer stellar disc region and is used to 
obtain the ellipticity and position angle of the outer stellar 
disc, edisciim and PAdisc u m respectively, from the ellipse free- 
fit e and PA radial profiles (see Fig. 2]). Fig. U shows the 
fitting procedure for an example galaxy. 

In some cases the fit failed at the stellar disc limit 
ddisc Um due to ellipse error limits in the radial intensity gra- 
dient. However, adjusting this threshold allowed these fits 
to be forced out to eidisciim ■ In order to ensure that the el- 
lipticity edisciim and position angle PAdisciim of the outer 
stellar disc are reliable, we compared the values obtained 
from the ellipse fit to a manual measurement of edisc Um and 
P4disciim measured directly off the ACS image during the 
visual inspection. The ellipse values agreed with the visually 
measured values in all cases. 
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However, in some cases the ellipse radial e and PA pro- 
files could be unstable at the stellar disc limit aaisciim- In 
these cases, 7 field (~ 5 per cent) and 13 cluster galaxies 
(~ 7 per cent), we use the last reliable e and PA obtained 
from the ellipse profiles as our estimate of e and PA for the 
outer stellar disc. 

In some other cases the free-fit failed completely at 
ddisciim- Subsequent inspection of the ACS image sometimes 
revealed clear structure at adisciim- In these cases, 27 field 
(~ 19 per cent) and 24 cluster galaxies (~ 13 per cent), we 
use the estimate of eaisciim and P-AcUsclim measured directly 
off the ACS image as our estimate of e and PA for the outer 
stellar disc. Unfortunately, in a few cases the ACS image 
was too noisy for a reliable manual measurement of eaisclim 
and P-Adisciim- In these cases, 7 field (~ 5 per cent) and 4 
cluster galaxies (~ 2 per cent), we use the value of the e 
and PA determined b y galfit in the STAGES Sersic mod- 
els (|Grav et a.1.1 l2009h as our estimate of e and PA for the 
outer stellar disc. 

A fixed-parameter fit (fixed centre, e, and PA) using the 
e and PA of the outer stellar disc obtained from the free-fit 
is then used to produce our final surface brightness profiles 
(see Fig. [4]). The fixed fits are forced to extend out well into 
the sky background (semi-major axis of 600 pixels). Dur- 
ing the fixed-fitting process we also perform 4 iterations of a 
3 — a rejection applied to deviant points below and above the 
average to smooth some of the bumps in the surface bright- 
ness profiles that are due to non-axisymmetric features, i.e. 
not part of the disc (e.g. star-forming regions and super- 
novae). The final profile fits then undergo a sky subtraction 
and surface brightness calibration in order to produce the 
final surface brightness profiles which will be discussed in 
the following sections. 

3.1 Photometric calibration 

The STAGES HST/ACS V-band images have pixel values 
in units of electrons n c . These can be converted to AB mag- 
nitudes per pixel using the following expression, 

/ PHOTFLAM x n c \ 
m AB = -2.51og 10 ( EX p TIME ) + ***** ^ 

where the AB zeropoint 

zptAn = PHOTZPT - 5log 10 (PHOTPLAM) + 18.6921. (3) 

The following FITS header keywords are used (ACS 
Datahandbook v5.fjf]), 

PHOTFLAM = 7.766405 x 10" 20 erg cm' 2 A" 1 electron" 1 , 
PHOTPLAM = 5.919369 x 10 3 A, 
and PHOTZPT = -21.10. 

The effective exposure time EXPTIME (after ACS sub- 
pixel drizzling) varies for each sample galaxy and is between 
700-770 seconds. The AB magnitude for each pixel mAB is 
then converted into a surface brightness /i using the ACS 
pixel scale (0.03"), 

fi — tuab + 2.5\og 10 (pixel scale 2 ). (4) 

2 Pavlovsky, C, et al. 2004, "ACS Data Handbook", Version 3.0, 
(Baltimore: STScI) 



This gives a surface brightness fi per pixel in 
magAB arcsec -2 in the V-band. 

Several corrections are also required in order to correct 
the surface brightness fi for galactic extinction, individual 
galaxy inclination i, and surface brightness dimming. 

(i) Galactic extinction: The STAGES field is affected by 
reddening d ue to substantial foreground dust (galactic ex- 
tinction, see lGrav et al.ll2009l) and an extinction correction 
^extinction of —0.18 magnitudes in the V^-band is required. 

(ii) Galaxy inclination: The surface brightness /i of a spi- 
ral galaxy increases with its inclination i to the line of sight. 
We correct for galaxy inclination using the surface bright- 
ness correction Ai = — 2.5 1og ln (l — e), wh ere e is the elliptic- 
ity of the outer isophote (|FreemarJll97(ih . We use the ellip- 
ticity determined for the fixed-fit edi sc m (outer isophote) to 
correct for galaxy inclination i. We do not attempt to correct 
the surface brightness for internal extinction (dust effects). 
However, since the galaxies studied here are all reasonably 
face-on ({q) — 0.7, (i) = 45°), and as we are dealing with the 
outer parts of the galactic disc, dust is expected to have very 
little effect on the surface brightness correction and no effect 
on the presence or strength of truncations/anti-truncations 
in the stellar disc. 

(iii) Surface brightness dimming: We correct the surface 
brightness /i profiles of our field galaxies so they are at 
the redshift of the cluster (z c \ = 0.167). This allows a fair 
comparison between galaxies across the redshift range of 
our sample. All cluster member galaxies are assumed to 
be at the cluster spectroscopic redshift. Surface brightness 
fl oc (1 + z)~ 4 for bolometric luminosity; however as we 
use the l/-band (F606W) filter (effective filter width depen- 
dent on z) we loose a (1 + z) term from this relation and 
(i oc (l + z)~ 3 . Therefore, correcting the /i of our field galax- 
ies requires multiplying the galaxy flux (in electrons n e ) by 
(1 + z ga i) 3 /(l + Zd) 3 , where z ga i is the photometric redshift 
of the field galaxy from COMBO- 17. This leads to a surface 
brightness correction A z for our field galaxies where 

^ = -2.51og 10 (|^). (5) 

Therefore, the final surface brightness /^corrected corrected 
for galactic extinction, galaxy inclination i, and surface 
brightness dimming is given by 

/^corrected = /^measured ^extinction 4™ Ai -(- A z . (6) 

3.2 Sky subtraction 

During the galfit Sersic model fitti ng performed by the 
GALAPAGOS pipeline (|Grav et alj|2009l ). the sky level is cal- 
culated individually for each source galaxy by evaluating a 
flux growth curve and using the full science frame. In this pa- 
per, for each sample galaxy we use the sky level determined 
by GALAPAGOS for our sky subtraction. Fig. [4] includes the 
free-fit radial intensity profile (in electron counts n e ) for an 
example galaxy showing the sky background determined by 
GALAPAGOS. 

The error in the GALAPAGOS sky level was estimated 
using observed 'dark' patches of sky located in some of the 
STAG ES ACS tiles (tiles 01, 10, 31, and 60, see lGrav et ail 
12009} ). Each patch (measuring 100 x 100 pixels) was visually 
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Figure 4. Example profile fit. Left column: ACS image (top) and bad pixel mask (bottom) showing the visually determined stellar disc 
limit a^isciim (black dashed line). The ACS image is shown with a logarithmic greyscale. Right column, top panel: The azimuthally 
averaged radial electron count n c profile from the free-parameter fit showing the GALAPAGOS sky background (red dashed line) and the 
stellar disc limit o-disc lim (black dashed line). Right column, second and third panels: The ellipticity e, and position angle PA radial 
profiles from the free-fit showing the position of the stellar disc limit a diac i im (black dashed line) used to obtain the e and PA of the 
outer stellar disc. PA is measured east of north. Right column, bottom panel: The final azimuthally averaged radial surface brightness 
profile from the fixed-parameter fit. The error in the surface brightness (grey shaded area) is due to over- and undersubtracting the sky 
background by ±1<t. The limiting surface brightness fj,\i m (blue dashed line) and critical surface brightness fi cv i t (blue solid line) represent 
the limit to which we trust profile breaks and the surface brightness profile respectively. 



Table 2. Determination of the ±lcr sky error from observed 
dark patches of sky in the STAGES ACS V-band imaging. 



STAGES tile 


-^Dark patches 


CT (n e ) (electron counts) 


01 


24 


0.1191 


10 


28 


0.1314 


31 


23 


0.3403 


60 


30 


0.1110 


±1<t sky error 


= Mean ct ( , 1c) 


0.1754 



inspected to ensure it was clear of any V-band sources. The 
distribution of electron counts n c (pixel values) in each dark 
patch is Gaussian in all cases. For each ACS tile we obtain 
the mean pixel value (n e ) in each dark patch of sky and 
then the standard deviation in the mean pixel values o~/ n<i \, 
see Table [2] The distribution of mean pixel values (n ) is 
approximately Gaussian for each tile. The final Tier sky 
error is the mean 07„ e \ from the 4 selected ACS tiles and is 
~ ±0.18 electron counts, see Tabled 

The sky subtraction error in the surface brightness pro- 
files due to the error in the GALAPAGOS sky background dom- 
inates over the individual errors produced by ellipse in the 



fitting process. At /i < 25 mag arcsec -2 the fit error domi- 
nates over the error in the sky subtraction but has a neg- 
ligible effect on the surface brightness profile. However, at 
jj. > 25 mag arcsec -2 the sky subtraction error dominates the 
error in the surface brightness profile. The sky subtraction 
error can have a significant effect on the surface brightness 
profile of the spiral galaxies, especially in the outer regions 
where the surface brightness fiv approaches that of the sky 
background. However, for any particular galaxy the global 
sky subtraction error is approximately constant across the 
length of the surface brightness profile. Therefore we can 
specify the error in our surface brightness profiles by gen- 
erating profiles for when the sky background is over- and 
undersubtracted by Tier (see Fig. [4j . 

The ±lcr error in the sky background corresponds to 
a critical surface brightness limit ^u, c j-tt below which the sky 
subtracted radial surface brightness profile of a galaxy be- 
comes unreliable. This critical surface brightness (Merit is 
approximately 27.7 mag arcsec -2 . We also define a limit- 
ing surface brightness /xn m , corresponding to a ±3cr sky er- 
ror, below which identifying profile breaks becomes unreli- 
able. The limiting surface brightness /in m is approximately 
26.5 mag arcsec 2 . 
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4 PROFILE CLASSIFICATION 
4.1 Profile inspection 

The azimuthally averaged radial surface brightness profile 
for each spiral galaxy in our field and cluster samples was 
visually inspected by three independent assessors (DTM, 
AAS, MEG) in order to identify potential profile breaks (in- 
flection points in the exponential region of the n profile). 
Three possible cases were considered: no break or simple ex- 
ponential profile, a single broken exponential, and cases with 
two profile breaks. In each case, break identification relates 
to the outer disc component of the galaxy radial surface 
brightness /j, profile and does not consider the inner varying 
bulge component. 

Our break identification is based solely on the surface 
brightness profiles and without direct inspection of the ACS 
images. We chose not to relate the profile breaks to visu- 
ally identified structural features because we wanted a break 
classification method that treated all galaxies equally, in a 
self-consistent manner, and avoided the prejudice that im- 
age inspection could introduce. Additionally, the aims of this 
work are to explore the effect of the environment on the 
structure of the outer regions of the galactic disc, regardless 
of the origins of any identified structural features. 

The radial limits of exponential regions either side of 
the break radius rb r k were also identified allowing for some 
bumps and wiggles due to spiral substructure (e.g. spiral 
arms and star forming regions). For the inner exponential 
disc, the inner boundary is chosen manually to avoid the 
region dominated by the bulge component. For the outer 
exponential disc, the outer boundary is generally taken to 
be where the surface brightness profile reaches the critical 
surface brightness (la above the sky background) but may 
be at higher \i depending on the nature of the profile. A small 
manually selected transition region (non-exponential) is al- 
lowed between exponential regions either side of the break. 
The break radius rbrk is defined as the mean radius for the 
radial limits of this transition region. 

Due to the subjective nature of some galaxy profile clas- 
sifications, the number of galaxies with either no, one, or two 
breaks varied subtly between the different assessors. To ac- 
count for this, in what follows we perform parallel analysis 
on the breaks identified by each assessor and compare the 
final results. 



4.2 Defining the Outer Stellar Disc 

In this study, we wish to consider the effect of the galaxy 
environment on the outer regions of the stellar disc for spiral 
galaxies. Therefore, we only wish to consider broken expo- 
nentials in the outer disc region of the [t, profile. Any effect 
of the environment should be stronger in these faint, frag- 
ile, outer regions. Fig. [S] shows the surface brightness at the 
break radius /Zbrk for galaxies with one, and two breaks. In 
the case of one break, for each assessor there appears to 
be a slight difference in the distribution of /Zbrk between 
the field and cluster samples with cluster breaks occurring 
at fainter surface brightnesses. Kolmogorov-Smirnov (K-S) 
tests between the respective field and cluster samples show 
that this environmental difference is driven by Type III (up- 
bending break, antitruncated) profiles with a significance at 



the 2-3<r level. We offer no interpretation of this result due to 
its uncertain nature. In the case of two inflection points, we 
consider the inner and outer breaks separately and the sepa- 
ration between the inner and outer break occurs at a surface 
brightness /i ~ 24 mag arcsec -2 (see Fig. [5|. The same con- 
clusion is reached for the break samples generated by each 
assessor. Inner and outer breaks may have different physi- 
cal origins. Therefore, in order to ensure we are comparing 
intrinsically similar fi breaks in the stellar disc, we limit 
the breaks analysed to those with /ibrk > 24 mag arcsec -2 
(outer disc breaks) and to the outermost break if two inflec- 
tion points are present in this range. There may be break 
features at surface brightnesses \i < 24 mag arcsec -2 , and 
these may be related to the galaxy environment. However, 
we focus on breaks in the outer disc as these are more likely 
to be susceptible to environmental effects. We acknowledge 
that in doing this there is some potential for missed envi- 
ronmental effects in the inner disc, therefore we do include 
some additional tests using the inner/initial break, see Sec- 
tion H31 

In Section 13.21 we defined a limiting surface brightness 
Miim (±3<7 sky error) of 26.5 mag arcsec -2 . In our analy- 
sis we do not trust breaks at fi fainter than the /^Hm level 
as this is where identification of breaks becomes unreliable 
due to the deviation of the profiles generated by over- and 
undersubtracting the sky by ±la (see Fig. |4}. We there- 
fore also restrict our analysis to profile breaks that have 
/Ubrk < 26.5 mag arcsec -2 . 

4.3 Profile Types 

We classify our galaxies into three main types, those clas- 
sified to be Type I , Type II , or Type III depending 
on break features in their outer stellar disc (o - outer, 
24 < /j,v < 26.5 mag arcsec -2 ) only. The classification used 
assumes only one break in the outer disc. If two breaks are 
present the outer break is used for classification. This is only 
the case for 4 (~ 3 per cent) field galaxies. If the profile has 
no break in the outer disc then the galaxy has a simple ex- 
ponential outer disc and is classified as Type I . Galaxies 
that are pure exponentials across the length of their surface 
brightness profile, aside from the varying bulge component 
(Type 1, a subset of Type I ), are also identified. If the 
profile has a down-break in their outer n profile then the 
galaxy has a truncation in the outer disc and is classified as 
Type II . However, if the profile has an up-break then the 
galaxy has an antitruncation in the outer disc and is clas- 
sified as Type III . Examples of each profile type (Type I, 
I , II , and III ) are shown in Fig. [6] along with their ACS 
images showing the break radius r^k isophote. Please note 
that this classificati on scheme is different t o that used by 
previous works fe.g. lPohlen fc Truiilloll200rj ). 

4.4 Determining scalelengths and break strength 

In the following, our exponential fits are obtained using a 
linear least-squares fit to the original profile between radii 
identified during the visual inspection (see Section [4. ip . 

For galaxies with no break in their outer disc (Type 
I /I, pure exponentials) we obtain the scalelength h using a 
simple exponential fit (h = 1.086 x Ar/A/ifi t (r)) across the 
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Figure 5. The distribution of break surface brightness ^tbrk- The surface brightness at the break radius ^tbrk f° r field (top row) and 
cluster (bottom row) galaxies as determined by DTM (left column), AAS (centre column), and MEG (right column). The distributions 
show galaxies with one break (grey shaded area), and both the inner (red line) and outer break (blue dashed line) of galaxies with two 
breaks. For the case of two breaks the separation between the inner and outer break occurs at a surface brightness of ~ 24magarcsec — 2 . 
The position of the break surface brightness cut 24 < /tb r k < 26.5 mag arcsec -2 used to create our final break samples is shown (black 
dashed lines). Respective sample sizes are shown in the legends. Contamination of the cluster sample by the field is < 25 per cent. 



length of the disc component. For bulgeless spiral galaxies 
with pure exponential discs (Type I Sd galaxies), the scale- 
length h was confirmed to be compatible with the GALFIT 
effective radius r e (r e ~ 1.7h). The mean random error in 
scalelength due to the exponential fitting routine is < 10 
per cent for these galaxies. The mean systematic error in 
the scalelength due to the error in the sky subtraction ±1<t 
(see Section [3. 2p is also < 10 per cent. 

For galaxies where an outer disc profile break was iden- 
tified (Type II /III ) , we obtain the scalelength h of expo- 
nential fits either side of the break radius rb r k, and therefore 
obtain an inner and outer exponential fit for the stellar disc. 
The inner exponential disc extends from a radius of 7"i n , m i n 
to and has a scalelength hi n given by 

i i no/ 1 7"in,max Tin min /—n 

h in = T086 x -, (7) 

Mfit(nn,max) — £ifit (rin,mm ) 

where /xg t is the surface brightness from the exponen- 
tial fit. Similarly, the outer exponential disc extends from a 
radius of r ou t,min to r ou t,max and has a scalelength h out given 
by 

/lout = 1.086 x r ° ut ' ma * ~ r ° ut ' mi " -. (8) 

J — Mflt (Tout . turn / 

For these galaxies, the mean random error in scalelength 
due to the exponential fitting routine is < 10 per cent for 
hi n and < 20 per cent for h ou t. In order to measure the 
strength of the profile break T we use the logarithm of the 
outer-to-inner scalelength ratio, 

T = log 10 ftout/fein. (9) 



The mean random error in T due to the exponential 
fitting routine is ~ ±0.1 (< 20 per cent). Fig. [7] shows the 
break strength T plotted against the outer scalelength /i ou t 
for Type II /IIIo galaxies (classified by DTM) with the er- 
rorbars in T due to the sky subtraction error ±lcr (see Sec- 
tion [3]2]). The mean systematic error in T due to the sky 
subtraction error is also ~ ±0.1. 

A Type I /I galaxy (pure exponential) has no break, 
and therefore has a break strength of T = 0. A Type II 
galaxy (down-bending break, truncation) has a smaller outer 
scalelength h ou t with respect to its inner scalelength hin, 
and therefore has a negative break strength (T < 0). Simi- 
larly, a Type III galaxy (up-bending break, antitruncation) 
has a larger outer scalelength h out with respect to its inner 
scalelength hi n , and therefore has a positive break strength 
(T > 0). We present a selection of our fi profiles with fitted 
exponential regions and overplotted break radii from both 
the field and cluster environments in Figs. [S] and [§] respec- 
tively. The ACS images and /i profiles for all galaxies in our 
samples are presented in Appendix A (online version only). 

In the outer regions of the surface brightness profile, 
negative sky-subtracted flux can occur. Surface brightness 
cannot be defined for a negative flux. Therefore, these points 
are removed from our linear exponential fits to surface 
brightness. However, this could potentially lead to a slight 
measurement bias and a larger outer scalelength h out for 
some profiles. In order to address this issue we repeated 
all measurements of scalelength (hi n and /i ou t) and break 
strength T using a non-linear least squares exponential fit 
to the flux radial profile (retaining negative fluxes). We then 
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Figure 6. Examples of each class of profile: Type I, Type I , Type II , and Type III (top to bottom). Left panels: ACS V-band 
images. Right panels: Azimuthally averaged V-band radial surface brightness profiles. We overplot the break radii in the outer stellar 
disc, 24 < fi < 26.5 mag arcsec -2 (red solid lines). Also shown are break radii with [ibrk < 24 mag arcsec -2 (red dashed lines), which 
are not considered in this study, and the stellar disc limit ddisclim (black dashed line). The inner and outer scalelength, h; n and h ou t 
respectively, and the break strength T are also shown for reference. The ACS images are in a logarithmic greyscale. 
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Figure 8. Example azimuthally averaged radial surface brightness profiles for different profile types in the field environment. First 
column: Type 1 profiles (pure exponential profiles). Second column: Type I profiles which have a simple exponential region in the outer 
disc (24 < fiy < 26.5 mag arcsec -2 , region between blue dashed and dotted lines). Third column: Type 1I profiles (down-bending break 
in outer disc, outer disc truncations). Forth column: Type III profiles (up-bending break in outer disc, outer disc antitruncations) . The 
error in the fi profiles is for an over- and undersubtraction of the sky by ±fo\ The critical surface brightness fi CI n (±fcr sky error, blue 
solid line) and limiting surface brightness /iii m (±3c sky error, blue dashed line) are shown on the profiles. The position of outer disc 
profile breaks (24 < /ibrk < 26.5 mag arcsec -2 ) are overplotted as red solid lines, breaks with /^brk < 24 mag arcsec 2 are represented by 
red dashed lines. The least square exponential fits to exponential regions are overplotted on the relevant sections of the surface brightness 
profiles. The Hubble type, redshift z, scalelength h, and break strength T of the galaxies are also shown for reference. The ACS images 
and fi profiles for all galaxies in our samples are presented in Appendix A (online version only). 
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Figure 9. Example azimuthally averaged radial surface brightness profiles for different profile types in the cluster environment. Figure 
same as Fig.|8]but for the cluster environment. The ACS images and fi profiles for all galaxies in our samples are presented in Appendix A 
(online version only). 
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Figure 7. A comparison of break strength T and outer scale- 
length hout for our Type II /III galaxies in the field (blue 
crosses) and cluster (red points) environment (as classified by 
DTM). An outlying cluster galaxy located at /i ou t = 453.2 kpc, 
T = 2.382 is not shown for clarity. The T errorbars represent the 
error from the sky subtraction (Her) on a galaxy-galaxy basis. 
The mean error in T is ±0.1. Respective sample sizes are shown 
in the legend. 

compared these results to those from the linear fits to sur- 
face brightness in order to assess if any measurement bias 
was present. The results from both fitting methods were in 
very good agreement and no bias was observed. Therefore, 
we use the linear fits to surface brightness in our analysis. 



5 RESULTS 

The frequency of profile classifications by each assessor 
(DTM, A AS, MEG) for spiral galaxies in the field and clus- 
ter environments are shown in Table [3] In each case the 
profile classification is based on \iy breaks in the outer stel- 
lar disc (24 < nv < 26.5 mag arcsec -2 ). Uncertainty in the 
frequency/fraction of profile types Sf (/ = Ni/Ntota\) are 
calculated using 

( 5 Ji V _ ( SN * \ 2 ( gjVtotaj \ 2 2V Nj TVtotal f , 

\fi ) \Ni) \ ATtotai ) AWtotai ' [ ' 

where SNi = y/Ni and (Wtotai = V-Wtotai . 

Due to the subjective nature of some profile classifi- 
cations, the frequency obtained for each profile type varies 
subtly between the different assessors. The agreement be- 
tween the three assessors is generally good for the outer disc 
profile types (Type I , II , and III ). This agreement is very 
good in the cluster, but slightly weaker in the field due to 
most of the galaxies being at higher redshift and therefore 
having poorer quality profiles. There is less agreement be- 
tween the three assessors for the frequency of Type I profiles 
(pure exponentials) due to the subjective nature of the clas- 



sification of weak profile breaks across the length of the n 
profile and inner profile breaks near the bulge component. 

The frequency of profile types in the outer stellar disc 
(Type I , II , and III ) are approximately the same in the 
field and cluster environment. For both field and cluster spi- 
rals, ~ 50 per cent have a simple exponential profile in the 
outer stellar disc (Type I ), ~ 10 per cent exhibit a down- 
bending break (outer disc truncation, Type II ), and ~ 40 
per cent exhibit an up-bending break (outer disc antitrun- 
cation, Type III ). The frequency of Type I profiles is also 
the same in the field and cluster for each assessor and is 
~ 20 ± 10 per cent. We stress that due to our profile classi- 
fication being based on breaks with fibrk > 24 mag arcsec -2 , 
our profile type fractions do not necessarily need to agre e 
with those of previous works (e.g. lPohlen fc Truiillol 12006). 
For our cluster galaxies, we also investigated the frequency 
of outer disc profile types as a function of s tellar surface mass 
density (range 10 10 -10 13 M Q Mpc" 2 , see IWolf et all (|2009h 
for details of how this was measured). However, no signif- 
icant differences were observed. These results suggest that 
the profile type in the outer disc of spiral galaxies is not 
affected by the environment from the field out to the inter- 
mediate densities of the A901/2 clusters. 

In this work, we define the outer disc to be 24 < /i < 
26.5 mag arcsec -2 , see Section 14.21 Our outer disc profile 
types are based on this somewhat arbitrary but justified 
definition. Independently and uniformly adjusting the upper 
and lower fi limits for the outer disc by +0.5 mag arcsec -2 
affects our profile type fractions by < 10 per cent. In each 
case, there remains no significant differences between the 
fractions in the field and cluster environment. 

The majority of our field sample have redshift z ~ 0.23 
while the cluster galaxies have z ~ 0.167. However, evolu- 
tionary effects are not expected to have an impact on the 
results of this work. The break strength T of our field galax- 
ies show no correlation with redshift and evolutionary effects 
on the scalelength h between the mean redshifts of our field 
and cluster samples i s only ~ 5 per cent (based on the fits of 
iBuitrago et alJfeOQSl for the expected size evolution of mas- 
sive disc galaxies). 

5.1 Morphology 

For each outer disc profile type (Type I , Ho, and III ), we 
compare the morphological mix (Sa, Sb, Sc, Sd) in the field 
and cluster environment to that of the total spiral galaxy 
sample (see Fig. I10|) . We also considered whether or not the 
galaxies were barred. 

The morphological mix of the total galaxy sample is ap- 
proximately the same in the field and cluster, with perhaps 
a slight excess of Sa galaxies in the cluster environment (Sa 
galaxies are visually defined to have a high bulge-to-disc 
(B/D) ratio, see Section [2]). This excess is dominated by 
barred Sa galaxies. For each outer disc profile type (Type 
Io, II , and III ), the morphological mix is also approxi- 
mately the same in the field and cluster environment. How- 
ever, there is a potential slight excess of Sa galaxies in the 
cluster for both Type II and III galaxies. In both cases 
the excess is dominated by barred Sa galaxies. Except for 
the slight excess of barred Sa galaxies in the cluster, the 
morphological mix for barred and unbarred galaxies are the 
same in the field and cluster for each profile type. 
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Table 3. The frequency of profile types in the field and cluster environments for the three independent assessors (DTM, AAS, MEG). 
Type I is a sub-sample of Type I Q and percentage errors are calculated using equation 1 101 



Assessor 


Type I 


Outer disc profile types 


Unclassified 






Type I 


Type II 


Type III 




Field galaxies 












DTM 


17 (12 ±2%) 


59 (41 ± 2 %) 


18 (12 ±2%) 


68 (47 ± 2 %) 





AAS 


42 (29 ± 3 %) 


76 (52 ± 2 %) 


14 (10 ±2%) 


55 (38 ± 2 %) 





MEG 


38 (26 ± 3 %) 


86 (59 ± 2 %) 


16 (11 ±2%) 


40 (28 ± 3 %) 


3 (2 ± 2 %) 


Cluster galaxies 












DTM 


28 (15 ± 2 %) 


80 (44 ± 2 %) 


19 (10 ±2%) 


83 (46 ± 2 %) 





AAS 


48 (26 ± 2 %) 


83 (46 ± 2 %) 


23 (13 ±2%) 


74 (41 ± 2 %) 


2 (1 ± 1 %) 


MEG 


48 (26 ± 2 %) 


88 (48 ± 2 %) 


17 (9 ±2%) 


69 (38 ± 2 %) 


8 (4 ± 2 %) 



The morphological mix of Type I Q galaxies is compa- 
rable to that of the total galaxy sample. This implies that 
Type I galaxies are equally probable for all morphological 
types. For Type II galaxies, the morphological mix frac- 
tions are affected by low number statistics. However, there 
is a clear trend with the Hubble type fraction (Sa-Sc) ris- 
ing towards later Hubble types. For Type III galaxies, the 
morphological mix is also comparable to that of the total 
galaxy sample. There is also a clear trend with the Hub- 
ble type fraction (Sa-Sd) decreasing towards later Hubble 
types. 

We conclude that outer disc truncations (Type II ) are 
more prevalent in later Hubble types while outer disc an- 
titruncations (Type III ) are more prevalent in early Hubble 
types. However, there is no evidence to suggest the morpho- 
logical mix for the different outer disc profile types is affected 
by the galaxy environment. 

We also compare the frequency of outer disc profile 
types (Type I , II , and III ) for different Hubble type mor- 
phologies (Sa, Sb, Sc, Sd) in both the field and cluster envi- 
ronment (see Fig. For Type I galaxies, we observe no 
clear correlation between the frequency of profile type and 
the Hubble type in both the field and cluster environment. 
The frequency of Type I galaxies is approximately the same 
for each Hubble type morphology, perhaps raising slightly 
towards later Hubble types in the cluster environment. For 
Type II and Type III galaxies, a weak trend is observed 
between the frequency of profile type and the Hubble type 
in both the field and cluster environments. The frequency 
of Type II galaxies increases slightly towards later Hubble 
types while the frequency of Type III galaxies decreases. 
For Sd galaxies, low number statistics could be masking the 
continuation of these trends. Th ese results are in qualita - 
tive agreement with the work of iPohlen k, Truiilld (|2006t ) 
who find a clear correlation between the frequency of profile 
type and Hubble type using a sample of ~ 90 nearby galax- 
ies from SDSS (classification based on fi breaks across the 
entire stellar disc so direct comparisons cannot be made). 
We conclude that there is no significant difference between 
the frequency of outer disc profile types for different Hubble 
type morphologies in the field and cluster environments. 

5.2 Pure exponential outer discs (Type I ) 

For galaxies where no outer disc fi break was identified 
(Type I galaxies, which include all Type I galaxies), we 



compare the outer scalelength h ou t distributions in the field 
and cluster environments to see if there is any evidence for 
an environmental dependence on outer scalelength h ou t (see 
Fig. ll2p . The outer scalelengths hout for these Type I galax- 
ies were transformed into intrinsic linear scales using a cos- 
mology of Ho = 70kms" 1 Mpc" 1 , fi A = 0.7, and ft m = 0.3. 
We use the fixed cluster redshift (z = 0.167) to determine 
the intrinsic outer scalelengths of our cluster Type I galax- 
ies and the original COMBO-17 redshift estimate for our 
field Type I galaxies. Therefore, photo-z errors only prop- 
agate into the intrinsic scalelengths of our field galaxies and 
not our cluster galaxies. The mean err or in h ou t associate d 
with the photo- z error is < 10 per cent (' Malt by et alj|20"icl ). 

In our surface brightness n profiles, the error in the 
GALAPAGOS sky (see Section [3~2)) can have a significant effect 
on the scalelength h and break strength T, especially in the 
outer regions where the surface brightness n approaches the 
critical surface brightness /i cr it (27.7 mag arcsec~ 2 ). How- 
ever, for any particular galaxy the sky subtraction error can 
be taken to be approximately constant across the length of 
the surface brightness profile. Therefore, we can account for 
this error by performing parallel analysis for when the sky 
background is over- and undersubtracted by ±lo\ The mean 
error in h ou t due to the sky subtraction error is ~ ±0.3 kpc 
(< 10 per cent) for our Type I galaxies. Random errors in 
h ut due to exponential fitting are also typically < 10 per 
cent (see Section f4.4p . We also perform parallel analysis on 
the Type I samples generated by the three assessors (DTM, 
AAS, MEG) in order to account for the subjective nature of 
the profile classifications and compare the final results. 

In all cases (all parallel analyses), we observe no clear 
difference between the distributions of outer disc scalelength 
h ou t for our Type I spiral galaxies in the field and cluster 
environments (see Fig. I12[) . 

In order to test the significance of these results we con- 
struct outer disc scalelength h ou t cumulative distribution 
functions (CDFs, see Fig |12|l for the Type I galaxy samples 
and perform Kolmogorov-Smirnov (K-S) tests between cor- 
responding samples from the field and cluster environments. 
These K-S tests are used in order to obtain the probability 
that the field and cluster Type I samples are not drawn 
from the same continuous h ou t distributions. The results of 
these K-S tests are shown in Table U 

In this work, we only consider an environmental effect 
on the Type I outer scalelength h ou t to be significant if K- 
S tests yield a 2er level probability for the field and cluster 
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Figure 12. Comparing outer disc scalelength h ou t distributions in different environments. Top row: Outer disc scalelength h ou t distribu- 
tions for Type 1 galaxies in the field (blue line) and cluster (red dashed line) environment as classified by DTM (left panel), AAS (centre 
panel), and MEG (right panel). Bottom row: The corresponding outer scalelength h ou t CDFs showing the probability p that compared 
samples are not drawn from the same continuous h ou t distributions in the bottom right of each plot. Respective sample sizes are shown 
in the legends. Random errors in scalelength are typically < 10 per cent. Systematic errors in scalelength due to the error in the sky 
subtraction are also typically < 10 per cent. Contamination of the cluster sample by the field is < 25 per cent. We find no significant 
difference between the CDFs in each environment and no evidence to suggest that the outer scalelength h ou t of our Type I galaxies are 
not drawn from the same continuous /i ou t distributions. 



Type I galaxy samples not being drawn from the same con- 
tinuous h out distributions. This probability p (Field/ Cluster) 
is below the 2a level for each assessor and for when the sky 
background is over- and undersubtracted by ±ler. Therefore, 
we find no evidence to suggest that the outer disc scalelength 
/tout of Type I galaxies is dependent on the galaxy environ- 
ment. 

We perform the same analysis for Type I galaxies, a 
sub-sample of Type I galaxies. These galaxies have simple 
exponential fj, profiles across the length of their disc com- 
ponent. The K-S test results for these galaxies are shown 
in Table |U The probability p(Field/Cluster) is below the 2a 
level in each case. Therefore, we also find no evidence to sug- 
gest that the scalelength h of Type I galaxies is dependent 
on the galaxy environment. 



5.3 Outer disc breaks (Type II /III ) 

For our Type II and III galaxies, we compare the distribu- 
tion of break strength T in the field and cluster environments 
to see if there is any evidence for an environmental depen- 
dence on break strength T in the outer disc (see Fig. I13[) . 
We perform similar parallel analyses and statistical tests as 
for our Type I galaxies. The mean error in T due to the 
sky subtraction error is ±0.1. Random errors in T due to 



exponential fitting are also typically ±0.1 (see Section H. 41) . 
The results of the K-S tests are shown in Table [4] 

In all cases (all assessor samples and sky versions), 
we observe no clear difference between the T distribu- 
tions in the field and cluster environments. The probability 
p(Field/Cluster) is below the 2a level in each case. There- 
fore, we find no evidence to suggest that the break strength 
T of our Type II /III galaxies is dependent on the galaxy 
environment and this result is robust to the error in the sky 
background and the subjective nature of the profile classifi- 
cations. 

We explore this result by splitting the combined field 
and cluster Type II /III galaxy sample into two sub- 
samples about the median value for stellar mass M* im d n 
(1O 1O M0). Similar statistical tests are then used to deter- 
mine the significance of an environmental dependence inde- 
pendently for high- (M» > M* >m dn) an d low-mass (M» < 
A^*,mdn) Type II /III galaxies. A similar test was also per- 
formed using B — V colour instead of stellar mass where B 
and V were obta ined from the STAGES master catalogue 
i|Grav et alj|2009h and (B- V), mdn = 0.70 mag. In all cases 
(all assessor samples and sky versions), the significance of 
an environmental dependence on break strength T remains 
below the 2a level. Therefore, we conclude that there is no 
evidence to suggest an environmental dependence on break 
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Figure 13. Comparing break strength T (log 10 /i ut/^i n ) distributions in different environments. Top row: Break strength T distributions 
for Type II /ni galaxies in the field {blue line) and cluster {red dashed line) environment as classified by DTM {left column), AAS 
{centre column), and MEG {right column). Bottom row: The corresponding break strength T CDFs showing the probability p that 
compared samples are not drawn from the same continuous T distributions in the bottom right of each plot. Respective sample sizes are 
shown in the legends. Random errors in break strength T are typically < 0.1. Systematic errors in break strength due to the error in 
the sky subtraction are also ~ ±0.1. Contamination of the cluster sample by the field is < 25 per cent. We find no significant difference 
between the CDFs in each environment and no evidence to suggest that the break strength T of our Type II /I1I galaxies are not drawn 
from the same continuous h ou t distributions. 



strength T is dependent on either the stellar mass or the 
B — V colour of Type II /IIIo galaxies. 

Another test performed was a comparison of the break 
radius 7- D rk (in units of the GALFIT effective radius) in the 
field and cluster environments for both our Type II (T < 0) 
and Type III (T > 0) galaxies separately. In both cases, no 
significant difference (< 2a) was observed between the break 
radius in the field and cluster environments. However, there 
is a slight indication that Type II breaks occur at smaller 
break radii than Type III breaks. 

We also compare the break strength T of our Type 
IIo/IIIo galaxies with the effective r adius r e determin ed by 
the STAGES galfit Sersic models l|Grav et al.l 120091 ) (see 
Fig. I14|l . The size-break-strength distribution is the same 
for both field and cluster galaxies. However, we note that 
at small effective radii (r e < r e , m dn) there is an absence of 
Type II galaxies (outer disc truncations, T < 0) in both the 
field and cluster environments. r e , m d n (3.58 kpc) is the me- 
dian effective radius of the combined field and cluster Type 
IIo/IIIo galaxy sample. In contrast, there is an abundance 
of Type III galaxies (outer disc antitruncations, T > 0) 
reaching down to r e ~ 1 kpc. A comparison of the break 
strength T distributions for large- (r e > r e , m< in) and small- 
r e (r e < r e , m d n ) galaxies shows a clear difference in the two 
distributions (see Fig. I14p . 

In order to test the significance of this result we con- 
struct break strength T CDFs (see Fig. [T4j , for large- and 



small-r e Type II /III galaxies and perform K-S tests be- 
tween the large- and small-r e galaxy samples to obtain the 
probability the different samples are not drawn from the 
same continuous break strength T distributions. Similar par- 
allel analyses are also performed as in previous tests. 

In all cases (all assessor samples and sky versions), 
the probability that the large- and small-r e galaxy sam- 
ples p(large-r e /small-r e ) are being drawn from different con- 
tinuous break strength T distributions is at the 3cr level. 
Therefore, we find some evidence to suggest that the break 
strength T of Type II /IIIo galaxies is dependent on the 
galaxy effective radius r e . We offer no interpretation of this 
result. 

5.4 Inner/initial disc breaks 

The focus of this paper is on the effect of the environment on 
the structure of the outer regions of galactic discs. By exam- 
ining only breaks in the outer disc {fi > 24 mag arcsec -2 ), 
we have tailored our analysis to allow for the assessment of 
an environmental effect on the outer regions of the galac- 
tic disc, where the environment is most likely to have an 
effect. However, we acknowledge that there is some poten- 
tial for missed environmental effects in breaks at a higher 
surface brightness (p, < 24 mag arcsec -2 ). Therefore, in or- 
der to test whether we could have missed any environmen- 
tal effect in the inner disc we repeated our analysis using 
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Figure 10. The morphological mix for outer disc profile types. 
The fraction of each Hubble type in the field (crosses) and cluster 
(circles) environment for (top to bottom) the total spiral galaxy 
sample, Type I galaxies, Type II galaxies, and Type III galax- 
ies. Respective sample sizes are shown in the legends and the er- 
rors in the Hubble Type fraction were calculated using equation 
1101 Contamination of the cluster sample by the field is < 25 per 
cent. We observe no significant difference between the morpho- 
logical mix in the field and cluster environment for each profile 
type. 



the inner/initial break for profile classification. However, we 
find no correlations of the frequency of profile type with the 
galaxy environment. We also find no significant effect of the 
environment on either the scalelength h of Type I galaxies 
(pure exponentials), or the break strength T of inner/initial 
breaks. 
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Figure 11. The fraction of outer disc profile types, Type I , 
(crosses), Type II (circles), and Type III (squares) for differ- 
ent spiral Hubble type morphologies in the field (top) and cluster 
(bottom) environment. The sample size for each profile type is 
shown in the legends and the sample size for each Hubble type is 
shown on the x axis. The errors in the profile type fraction were 
calculated using equation 1101 Contamination of the cluster sam- 
ple by the field is < 25 per cent. We observe similar trends for 
both field and cluster galaxies. The fraction of Type I is approxi- 
mately constant with Hubble type, perhaps rising slightly towards 
later Hubble types in the cluster. However, the fraction of Type 
Ho galaxies (truncations) increases slightly towards later Hubble 
types while the fraction of Type III galaxies (antitruncations) 
decreases. 



6 CONCLUSIONS 

We present an analysis of l/-band radial surface brightness 
profiles for spiral galaxies from the field and cluster envi- 
ronments using HST/ ACS imaging and data from STAGES. 
Using a large, mass-limited (M» > 1O 9 M0), visually clas- 
sified sample of ~ 330 field and cluster spiral galaxies, 
we assess the effect of the galaxy environment on the ra- 
dial surface brightness fi profile in the outer stellar disc 
(jj, > 24magarcsec~ 2 ). 

We classify our spiral galaxies according to fi break fea- 
tures in their outer disc. If the galaxy has no break in this [i 
range then the galaxy has a simple exponential profile in the 
outer disc and is classified as Type I . Galaxies that have 
a simple exponential profile across the length of their disc 
component (Type I, a subset of Type I ), are also identi- 
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Figure 14. A comparison between the break strength T 
(log 10 h ou t/hi n ) and physical size r e of our Type II /III galax- 
ies (denned by DTM). Top Panel: The effective radius re-break- 
strength T distributions for our field {blue crosses) and cluster 
(red points) Type II /III galaxies, showing the position of the 
median effective radius r e m( j n (3.58 kpc) of the combined field 
and cluster sample (black dashed line). An outlying cluster galaxy, 
located at T = 2.382, r e = 4.442 kpc is not shown for clarity. Mid- 
dle panel: A comparison of the distribution of break strength T 
for small-r e , r e < r e m( j n (blue line) and large-r e , r e > r e m< j n 
(red dashed line) Type II /IIIo galaxies. Bottom panel: A com- 
parison of the break strength T CDFs for small-r e (blue line) and 
large-r e (red dashed line) Type II /III galaxies. The probabil- 
ity p that small-r e and large-r e Type II /III samples are not 
drawn from the same continuous break strength T distributions 
is shown in the top left of the plot. Respective sample sizes are 
shown in the legend. We find a significant difference between the 
CDFs of the small-r e and large-r e Type II /III samples with a 
high probability (at the 3<r level) that the samples are not drawn 
from the same continuous break strength T distributions. 



Table 4. The K-S test results for Type I, Type I D , and Type 
II /III galaxies as classified by DTM, AAS, and MEG. K-S 
tests give the probability that the respective field and cluster 
p(Field/Cluster) samples are not drawn from the same continuous 
h distributions for Type I and Type I galaxies and T distribu- 
tions for Type II /IIIo galaxies. Results are also shown for when 
the sky is over- and undersubtracted by ±lcr. We find no environ- 
mental dependence on either the scalelength h of Type I and Type 
Io galaxies or the break strength T of Type II /III galaxies. 



p(Field/Cluster) 



Sky subtraction 


Under (-la) 


Nominal 


Over (+ 


Type I (h) 








DTM 


0.078 


0.101 


0.167 


AAS 


0.825 


0.761 


0.865 


MEG 


0.952 


0.919 


0.921 


Type Io (h) 








DTM 


0.096 


0.146 


0.186 


AAS 


0.545 


0.545 


0.859 


MEG 


0.720 


0.716 


0.629 


Type IIo/IIIo (T) 








DTM 


0.165 


0.038 


0.122 


AAS 


0.797 


0.908 


0.795 


MEG 


0.180 


0.307 


0.563 



la) 



fied. However, if the galaxy exhibits a broken exponential in 
the outer disc it is classified as either Type II for a down- 
bending break (outer disc truncation) or Type III for an 
up-bending break (outer disc antitruncation). 

The frequency of outer disc profile types (Type I , II , 
and III ) is approximately the same in both the field and 
cluster environment. For both field and cluster spirals, ~ 50 
per cent have a simple exponential profile in the outer stellar 
disc (Type I ), ~ 10 per cent exhibit a down-bending break 
(outer disc truncation, Type Ho), and ~ 40 per cent ex- 
hibit an up-bending break (outer disc antitruncation, Type 
IIIo). These results imply that the shape of the outer disc 
surface brightness profile is not dependent on the galaxy 
environment. These results agree for break classifications 
performed by three independent assessors. We stress that 
due to our profile classification being based on breaks with 
^ibrk > 24 mag arcsec -2 , our profile type fractions do not 
necessarily need to agre e with those of previous works (e.g. 
iPohlen fc Truiilldl2006l ). 

We also find that the morphological mix (Sa, Sb, Sc, 
Sd) of the different outer disc profile types is the approxi- 
mately same in the field and cluster environment. However, 
we do find a dependence of the shape of the fi profile in the 
outer stellar disc on the Hubble type. Outer disc truncations 
are slightly more frequent in later Hubble types while the 
outer disc antitruncations are slightly more frequent in ear- 
lier Hubble types. The same dependence is observed in both 
the field and the cluster envir onments. This result is in qual- 
itative agreement with that of lPohlen fc Tru ullo ( 2006j) who 
find a similar correlation using galaxies from SDSS (classifi- 
cation based on entire disc component so direct comparisons 
cannot be made). 

For galaxies with a purely exponential outer stellar disc 
(Type I , ~ 50 per cent), we find no evidence to suggest 
an environmental dependence on the outer disc scalelength 
ftout- We also find no evidence for an environmental depen- 
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dence on the scalelength h for galaxies that have a pure 
exponential profile across the length of their disc compo- 
nent (Type I, a subset of Type I ). For galaxies with a 
broken exponential in their outer stellar disc, either down- 
bending (Type IIo, outer disc truncation, ~ 10 per cent) 
or up-bending (Type III , outer disc antitruncation, ~ 40 
per cent), we measure the break strength as the outer-to- 
inner scalelength ratio log 10 h ut/hi n - We also find no evi- 
dence to suggest an environmental dependence on this break 
strength. We conclude that there is no evidence to suggest 
the stellar distribution in the outer stellar disc is affected 
by the galaxy environment for these galaxies. These results 
have been shown to be robust to the sky subtraction error 
and to agree for break classifications performed by three in- 
dependent assess ors. This work is in quali tative agreement 
with the work of IPohlen fc Truiillol |2006) who come to a 
similar conclusion for a sample of ~ 90 spiral galaxies from 
SDSS (classification based on entire disc component). 

We also find that for galaxies with small effective radii 
(r e < 3 kpc) , Type Il profiles (outer disc truncations) are 
rare in both the field and cluster environments. In contrast 
Type III (outer disc antitruncations) are commonplace. 

Our results suggest that the galaxy environment has no 
effect on the stellar distribution in the outer stellar disc of 
spiral galaxies from the field out to the intermediate densi- 
ties of the A901/2 clusters. This implies that the origin of 
broken exponentials is related to an internal mechanism (e.g. 
star formation threshold, resonance phenomenon) or minor 
mergers. Our results are consistent with previous work car- 
ried out on the effect of the g alaxy environm e nt on disc 
features in the STAGES survey. iMarinova et all (|2009h find 
that the optical fraction of bars among disc galaxies does 
not show evidence for any strong variation between the field 
and the A901/2 clusters, suggesting that the mass redistri- 
bution associated with bar formation within galactic discs 
is not a strong function of environment from the field out to 
intermediate densities. However, it is important to note that 
STAGES only covers an intermediate density environment 
(projected galaxy num ber density up to ~ 1600 gal Mpc -3 , 
iHeiderman et af1l2009l ) and not a high densit y environment 
(e.g the COMA cluster, ~ 10 4 gal Mpc -3 . iThe fc White! 
Il986f l. Therefore, it is important to investigate the effect 
of the environment on the outer disc structure at higher 
densities (e.g. the COMA cluster). 

The results presented here are for one survey field and 
one multicluster complex: therefore it is also important to 
investigate the environmental dependence of the outer disc 
structure of spiral galaxies in different surveys and across a 
wide redshift range. The results of such studies will enable 
the evolution of the outer disc to be investigated in different 
galaxy environments. Also the investigation into the origin 
of broken exponentials will depend on multiwavelength data 
being used to create colour profiles and stellar mass distri- 
butions. This will allow an assessment of whether the origin 
of broken exponentials is due to the distribution of stellar 
mass or a radial change in the age of the stellar population. 
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